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ABSTRACT 

We discuss the dynamics of dust grains subjected to uncompensated torques 
arising from H2 formation. In particular, we discuss grain dynamics when a 
grain spins down and goes through a "crossover" . As first pointed out by Spitzer 
& McGlynn (1979), the grain angular momentum before and after a crossover 
event are correlated, and the degree of this correlation affects the alignment of 
dust grains by paramagnetic dissipation. We calculate the correlation including 
the important effects of thermal fluctuations within the grain material. These 
fluctuations limit the degree to which the grain angular momentum J is coupled 
with the grain principal axis ai of maximal inertia. We show that this imperfect 
coupling of ai with J plays a critical role during crossovers and can substantially 
increase the efficiency of paramagnetic alignment for grains larger than 10~^ cm. 
As a result, we show that for reasonable choices of parameters, the observed 
alignment of a 10^^ cm grains could be produced by paramagnetic dissipation 
in suprathermally rotating grains, if radiative torques due to starlight were not 
present. We also show that the efficiency of mechanical alignment in the limit 
of long alignment times is not altered by the thermal fluctuations in the grain 
material. 



Subject headings: ISM: Magnetic field. Dust, Extinction - Polarization 



1. Introduction 

Understanding of the observed alignment of interstellar grains is a challenging problem 
of nearly a half century's standing (see Roberge 1996). Lacking a proper understanding of 
the alignment processes, we can only tentatively interpret polarimetric data in terms of the 
magnetic field. Indeed, polarizing grains can be aligned with long axes either perpendicular 
or parallel to the magnetic field, depending on what causes the alignment (see Lazarian 
1994); they can also be not aligned at all (see Goodman 1996). 
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One of the essential features of grain dynamics in the diffuse interstellar medium 
(henceforth ISM) is suprathermal rotation (Purcell 1975, 1979). Originally, three separate 
causes of suprathermal rotation were suggested: inelastic scattering of impinging atoms 
when the gas and grain temperature differ, photoelectric emission, and H2 formation on 
grain surfaces. The latter was shown to dominate the other two for typical conditions in 
the diffuse ISM (Purcell 1979). More recently, radiative torques due to starlight have been 
identified as a major mechanism driving suprathermal rotation (Draine & Weingartner 
1996, 1997). 

Alignment of grains rotating suprathermally differs considerably from the alignment of 
thermally rotating grains.]^ The theory of paramagnetic alignment of suprathermal grains 
was discussed by Purcell (1979) and Spitzer & McGlynn (1979); henceforth SM), and has 
been elaborated by Lazarian (1995b, c, 1996b). Until recently, mechanical alignment, e.g. 
alignment caused by a gaseous flow, was thought not to be applicable to suprathermally 
rotating grains, as rapid rotation makes the grains not susceptible to such a process. 
However, two new mechanisms of mechanical alignment of suprathermally rotating grains, 
namely, the "crossover" and "cross-section" mechanisms, have been suggested recently 
(Lazarian 1995d) and shown to be effective in interstellar regions with gas-grain streaming 
(Lazarian & Efroimsky 1996; Lazarian et al. 1996). 

The crossover event is the most important period in the dynamics of suprathermally 
rotating grains. The H2 formation sites on a grain surface have a finite "lifetime" t^, 
which may be determined by the "resurfacing" of the grain by accreted atoms (SM, Purcell 
1979) or poisoning of active sites by oxygen (Lazarian 1995c). Because of the changes in 
the resulting torque, the spin-up has a finite duration and this limits the paramagnetic 
alignment attainable. In the case of mechanical alignment, it is during the crossovers that 
the grain is susceptible to alignment due to gaseous bombardment. 

The pioneering study of SM showed that the direction of angular momentum before 
and after crossover are correlated. However, SM found that the correlation was insufficient 
for the paramagnetic mechanism to achieve significant alignment within the model they 
adopted. 

Recent progress in the understanding of certain subtle issues of grain dynamics has 
led us to reexamine the crossover process. SM assumed that during suprathermal rotation 



^See Davis & Greenstein (1951), Jones & Spitzer (1967), Mathis (1986), Roberge et al. (1993), Lazarian 
(1995a) for paramagnetic alignment of thermally rotating grains; and Gold (1951), Lazarian (1994), Roberge 
et al. (1995), Lazarian (1996a) for mechanical alignment of thermally rotating grains. 
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the grain angular momentum is perfectly aligned with the axis of major inertia^. This 
coupling arises from internal dissipation^ which, as known from theoretical mechanics, 
causes a spinning solid body to rotate about its axis of major inertia, this being the state 
of minimum rotational kinetic energy for fixed angular momentum. The assumption of 
perfect relaxation seems natural, as the time-scale for internal relaxation for suprathermally 
rotating grains is many orders of magnitude less than the time-scale of the spin-up, but it 
is not exact: it disregards thermal fluctuations within the grain body. 

In fact, it was shown in Lazarian (1994) that due to thermal fluctuations, the coupling 
above is never perfect (the formal theory of this phenomenon is elaborated in Lazarian & 
Roberge 1996). The component of angular momentum perpendicular to the axis of the 
major inertia, although tiny compared to the grain angular momentum during suprathermal 
rotation, is very important in the course of a crossover. We therefore reconsider the SM 
theory of crossovers in order to allow for the effects of thermal fluctuations. 

In §2 we pose the problem and present the necessary facts concerning incomplete 
internal relaxation. In §3 we derive the disorientation factor accounting both for thermal 
fluctuations within the grain material and for the effects of gaseous bombardment. The 
latter effect is of secondary importance in diffuse clouds but may be important in molecular 
clouds. The consequences of the incomplete disorientation on paramagnetic and mechanical 
alignment are discussed in §4 and the conclusions are presented in §5. 

2. The problem 

A crossover is the event that occurs between two sequential spin-ups when the 
component of J parallel to the axis of major inertia passes through zero. This is a critical 
period for grain dynamics, and during the crossover the grain is susceptible to disorientation, 
which will limit the effectiveness of paramagnetic alignment. If the grain is situated in a 
region with substantial gas-grain streaming, the grain is susceptible to mechanical alignment 
during crossovers. Our task in the present paper is to describe the evolution of grains 
through crossover events, accounting for the effects of thermal fluctuations within the grain 
material. 

To understand the crossover one needs to recall certain basic features of suprathermal 
rotation. Here we assume that grains are spun-up by torques arising from H2 formation, and 



^For brevity, we refer to the principal axis of largest moment of inertia as the "axis of major inertia". 
■^The most important internal dissipation process is Barnett relaxation (Purcell 1979). 
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consider a "brick" with dimensions b x b x a and density p^. The ratio r = b/2a determines 
the degree of grain oblateness; r = 1 for the grain 2:2:1 discussed in Purcell (1979). 
It is possible to show (SM; Draine & Lazarian 1996) that the components of the torque 
perpendicular to the axis of major inertia average out and therefore only the component of 
the torque parallel to this axis matters. We direct the z-axis along the axis of major inertia. 
We let 

nn = n{}l)+2n{R2) , (1) 

where n(H) and n{B.2) are the concentrations of atomic and molecular hydrogen, 
respectively; the H2 fraction we denote y = 2n(H2)/nH. 

The number of H2 molecules ejected per second from an individual site is 
7a^nHfH(l — y)r{r + l)i^~^, where 7 is the fraction of H atoms (with mean speed and 
mass mn) adsorbed by the grain, and u is the number of active sites over the grain surface. 
The mean square torque from H2 formation is (see Appendix A) 

(Ll) ^ ll'a%l - yfnlmnvlEv-\\r + 1) , (2) 

where E ^ 0.2 eV is the kinetic energy of a nascent H2 molecule. The fluctuating torque Lz 
spins up grains to an rms angular velocity 

(Purcell 1979), where Iz = Ip^a^r^ is the z component of the momentum of inertia, td is the 
rotational damping time (see Appendix A) 

(r + 2) nHmHi'H (1-2 - 0.292?/) ' ^ ^ 

and is the lifetime of an active site. 

To obtain both characteristic numerical values and functional dependencies we 
will use quantities normalized by their standard values (see Table 1). We denote the 
normalized values by symbols with hats, e.g. d = a/(10~^ cm), with 10~^ cm as the 
standard value of grain size. We consider 'standard' an H2 formation efficiency 7 = 0.2 and 
fH = 1-5 ■ 10^ cm s~^. For diffuse clouds we assume that all hydrogen there is in atomic 
form and therefore y = 0. Note that sometimes the choice of 'standard' values is somewhat 
arbitrary, e.g. for the time being we assume the density of active sitesQ to be 10^^ cm~^, 
so that z/ = 80am(i^r(r + 1). 



"^The density of active sites depends on the interplay of the processes of photodesorption and poisoning. 
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Using standard values of the parameters we obtain the following expression for the 
angular velocity 



VamH/ 



\amnJ lo? (1.2 - 0.2927/) r3/2(r + 2) Vtd + 



Q;V2a2 2 - 0.2927/) r^l'^{r + 2) Vtd + II 

The lifetime of an active site of H2 formation is limited by both accretion of a mono-layer of 
refractory material (SM) and poisoning by atomic oxygen (Lazarian 1995c). The former is 
usually slower than the latter and could provide long-lived spin up with ti, 3> t^. Further on 
we use the term "resurfacing" to refer to the fastest mechanism of the two. The component 
of the mean torque along the axis of major inertia before and after, say, resurfacing may be 
directed either in the same direction as before the process, or in the opposite direction. In 
the latter case the grain undergoes a spin-down. 

The mean interval between crossovers is (Purcell 1979) 

t,^T^{tLtdf'^ . (6) 

The active site lifetime II is very uncertain. In our numerical examples in this paper we 
will take II — 10^^ s. For our standard parameters (Table 2) in a diffuse HI cloud this 
corresponds to tL/ta — 0.25/a, and ~ 1.6a^/^td. 

When a grain rotates about an arbitrary axis, the angular velocity precesses in grain 
body coordinates. The Barnett effect produces a magnetic moment parallel to the angular 
velocity. Purcell was the first to realize that this should result in internal dissipation with a 
dissipation time-scale inversely proportional to the angular velocity squared (Purcell 1979). 
It is possible to show (see §3.2) that this effect suppresses rotation around any axis but the 
axis of major inertia on a time-scale Ib ~ 10^/K s, where K is the ratio of grain rotational 
energy to the equipartition energy ~ kT . Both H2 formation (Purcell 1979) and radiative 
torques (Draine & Weingartncr 1996a) can produce K > 10^ and therefore grains tend to 
rotate around their axes of major inertia. 



It is shown in Lazarian (1995c, 1996b) that the poisoning intensifies when the number of active sites becomes 
greater than a critical number. The latter number depends on the migration time-scale of hydrogen and 
the activation barrier for the reaction between physically and chemically adsorbed hydrogen atoms. As 
the details of the grain chemistry are poorly known the critical number of active sites is highly uncertain. 
Therefore, for the sake of simplicity, we do not discuss differences in poisoning of active sites whenever their 
number exceeds the critical number and assume, following SM, a constant surface density a of active sites. 
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Although the ratio t^/tx ~ tB/td can be as small as 10~^ (see Table 2) the alignment 
of J with the axis of major inertia (ai) is not perfect. The deviations of ai from J arise 
from thermal fluctuations within the grain material (Lazarian 1994, Lazarian & Roberge 
1996). To estimate the value of such deviations recall that rotation about ai corresponds 
to the minimum of the grain kinetic energy for fixed J (internal dissipation does not alter 
J)0. For a symmetric oblate grain with 1^ > Ix = ly = I± (i-e., our "brick" with b > a), the 
grain kinetic energy is 

Ek{P) = ^(l + sm'p(^-l]) , (7) 



2/. V V/_ 

where (3 is the angle between J and ai. In thermodynamic equilibrium the fluctuations of 
the kinetic energy should have a Boltzmann distribution: 

/(/?)d/? = const -sin/Jexpf^-^^^d/J , (8) 

where is the dust temperature. It follows from (^ that the fluctuating component 
of angular momentum perpendicular to the axis of the major inertia {J]_) <^ can be 
approximated 

We may define the "thermal transverse angular velocity" 

As we will see below, lo\s) is the characteristic value for the minimum value of the grain 
angular velocity during a crossover. 

When the rotation is suprathermal, is negligible compared to J and angle /3 is 

very close to zero. However, as the component of angular momentum Jy parallel to the axis 
of major inertia decreases during crossovers, the angle /? = arctan( Jj_/ Jy) increases. 

When the angular velocity decreases sufficiently, internal relaxation becomes less 
efficient and the value of ( J^) rises as a result of the stochastic character of H2 formation 
and impacting gas atoms. At some point the component of J parallel to z passes through 
zero and the grain flips over (SM). Our task is to calculate the correlation of the grain 
angular momentum before and after the crossover. This is done in the next section. 



^If J does change due to H2 torques one should substitute its value averaged over the time of internal 
relaxation in Eq.(|^). 
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3. Crossovers 

In our treatment below we repeat the reasoning introduced in SM but with allowance 
for thermal fluctuations. The zeroth approximation, following SM, is the dynamics of a 
grain subjected to regular torques only. The dynamical effects of the stochastic torques can 
be evaluated by an approximate theory based on small perturbations of the zeroth-order 
solution. 



3.1. Zeroth approximation 

Let xyz be a coordinate system frozen into the grain, with z along the grain axis of 
major inertia ai. Let xoyoZo be an inertial coordinate system, with zo|| J (at some initial 
time). Let (3 be the angle between the z-axis and J: Jz = J cos (3. If no external torques 
act, then J = constant and the z-axis and the angular velocity w will each precess around J 
(or zo) at a frequency Up = (/^ - I^)JjIJz. 

Now consider the effect of a (weak) torque L which is fixed in body coordinates xyz. 
On time scales long compared to the rotation of the grain around w and the precession 
of w around zq imply that the only torque component which does not average to zero is 
that due to Lz, the component of L along the z-axis. After this averaging we obtain 

^ = L,cos/?i^. (11) 
From the Euler equations (see SM) we find the components of w in body coordinates 

= - = (12) 

u± = Jj_/ 1± = const, (13) 
where t = at the moment of crossover. Eq. (|TT]) can be rewritten 

-L,cos/5 = L,^ . (14) 



dt Jzo 



Since J. = Ld we obtain 



Ji=Lit' + Ji . (15) 

According to (|Il[) , the direction of J does not change - the torque L acts only to change 
its magnitude (see Eq. (p!5D). Thus the zeroth approximation predicts a perfect correlation 
between the angular momentum directions prior to and after the crossover. The stochastic 
torques make the story more involved. 
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3.2. Crossovers & Barnett fluctuations 

Our considerations above ignored the fluctuations associated with the Barnett effect. 
As discussed by Lazarian & Roberge (1996), angle /? fluctuates on the time of the Barnett 
relaxation (Purcell 1979): 

tB{uj) = ^ , (16) 

where 

^ = 7.1 X 10i^a2^frfi^o(r) . (17) 

and 

3 [Ar^ + 1)3 

*°W-125;W^ ■ '''' 
Note that the ratio r = 1/2 corresponds to a cubical grain, for which no internal relaxation 
is expected in agreement with Eq. (117 



The fluctuations in (3 span the interval (0, tt) when J J±q (Lazarian & Roberge 
1996). SM showed that during crossovers J J± and therefore such fluctuations must be 
accounted for provided that tc > tB{uj±), where the crossover time is 

tc=^ , (19) 

where J is the time derivative of J. 

When ^^(ti;^) ^ Barnett fluctuations will cause /5 to range over the interval (0,7r), 
resulting in frequent reversals of the torque in inertial coordinates. Consequently, the actual 
time spent during the crossover will be increased. Quantitative analysis of this regime is 
beyond the scope of the present paper; it does appear clear, however, that each crossover 
will be accompanied by substantial disalignment when ts '^tc- 

In our present study we conflne ourselves to the other limiting case, namely, tc/ts <^ 1, 
in which case the Barnett fluctuations during a crossover can be disregarded and the 
initial distribution of J± with the mean value given by Eq. (^ is produced by the Barnett 
fluctuations during the long time interval between crossovers. 

We shall prove below that for typical interstellar conditions the value of J± mostly 
arises from thermal fluctuations within the grain material during intervals of suprathermal 
rotation and therefore is given by Eq. (^). Thus we can estimate tc- 

t^^?m^2.6.10'(^f \. M r)^ . (20) 

J \ E J -fUnVHil - y) 



9 



where 




3r(l + 4r2) 



5(4r2 - 1) 




(21) 



The ratio 



tB 

tc 

where Oc is the critical radius Oc, which 




(22) 



1.47 X 10 



-5 



( 



(1 - yfn^jv'jjEpl 




cm 



(23) 



where 




234375r5 



(24) 



(4r2 



l)(4r2 + l)7 



It follows from our discussion above that we attempt to deal only with the case a > Oc 
while leaving the more complex regime a < Qc, where Barnett fluctuations during the 
crossover are important, to be dealt with elsewhere. We use the inequality a > ac rather 
than a ^ ttc due to the strong dependence of the time ratio tc/ts on a: for a = 10"^ cm 
tc > 12tB, while for a = 2.0 x 10~^ cm ts > 7tc- Below we analyze the implications of the 
critical size Oc in the context of the variations of alignment with grain size. 

The numerical value ac ~ 1.5 x 10^^ cm is quite robust: the most uncertain grain 
parameter is the surface density of active sites a, but even varying a by a factor 10^ changes 
Oc by only a factor 10^/^^ ~ 1.36. 

Although so far we have considered only suprathermal rotation driven by H2 formation, 
the existence of a critical size seems to be generic to the problem of disorientation in the 
course of crossovers. 

In our treatment above we disregarded gaseous friction. This is a good approximation 
in the diffuse medium since tc -C td (see Table 2). In molecular clouds as ?/ — 1 the two 
time scales may become comparable (e.g. if the atomic hydrogen fraction (1 — ?/) drops 
below 10^^), and gas drag should be included. We also assume t^ ^ tc, in the case of 
ti "^tc it becomes important to allow for variations in the time-averaged torque during 
the crossover. 



3.3. Random torques 



We consider the dynamical evolution given by Eq. (|T5|) as a zeroth-order solution of 
the problem, and the dynamical effects produced by stochastic torques as perturbations of 
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this solution. 

Each torque event produces an impulsive change of angular momentum 
A J = A JjZ + AJj^. The angular deviations of J in the (J, z)-plane are given 

by 

A r;|| = -^J^^^^f^ + ^J^^^o^^ (25) 



and in the transverse direction by 



A,. = ^ , (26) 



where AJ^i and AJ±2 are the components of AJ^ in the plane parallel and perpendicular 
to the (J, z)-plane, respectively. The grain is subject to the action of various torques. 
Here we discuss only torques arising from H2 formation (A J)h and gaseous bombardment, 
{AJ)g. Let Ni be the rate of H2 formation, and A^2 be the rate of gas-grain collisions. To 
simplify our notation, we denote the mean square change in angular momentum per H2 
formation event 

(AJ,)2 = (AJ,)| + ^(AJ,)J , (27) 

where N2/N1 is the number of gas-grain collisions per H2 formation event. The mean 
quadratic deviation of J per H2 formation is 

{A{rjf) = J-^{((AJ,)2)sin2/3+((AJxi)^)cos^/?+((AJ^2)')-(AJ,AJxi)sin2/?} , (28) 

where (A A J±i) — due to rotation around the axis of major inertia. For a symmetric 
oblate grain 

{{AJ^^f) = ((AJ^2)^) = liiAJ^f) , (29) 

and we obtain 

((Ar,)^) = J-|((AJ.)^)sin^/?+((AJ^)^)(l±|^)} . (30) 

The cumulative deflection due to i -|- 1 random impulses may be expressed as 

cos rji+i — cos rji cos( Ar^j) -|- sin rji sin( Ar^j) cos (pi , (31) 

where (pi is the angle between the deviation {Ar]i) and the great circle measured by rji. 
Averaging provides us with the result (see SM) 

(cos77/) = nj(cos A77i) . (32) 
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As Arj is small, we can expand cos Arj to obtain 



Hi exp 



Hence 



{cosrjf) ^ exp(— F) 
where F is the disorientation parameter 



F=- 
2 



N,{iArif)dt 



(33) 
(34) 

(35) 
and (|30|) 



and A^i = Lz/ (A J^) is the number of H2 torque events per second. From Eq. ( 
it is seen that both H2 torques and gaseous bombardment are included in Eq. 

To evaluate F, we obtain (3{t) from the zeroth-order grain dynamics with only regular 
torques {J± =const, = Lzt)] from tan/5 = J^/Jz we obtain 



dt 



dJ. 



Jxd/5 



J2 



-d/3 



Lz L^sin (3 J±Lz 
Substituting this into Eq. (|35| ) and integrating from /3 = to /5 = vr, one obtains 

^ TT {{AJzf) / , 3{{AJ^r)\ 



(36) 



Using Eq. (^TD and 



^J±\{AJz 
WM, we find 



1 



2 ((AJ,)2) 



(37) 



^^179. ((A J.);) 



224 Jx|(AJ^ 



where 



i^2(r) 



56 
179 



1 + 



3[16r2(r + 1) + 6r + 3] 
8r2(2r + 5) 



(38) 



(39) 



Eq. ( 137|) was derived assuming J_l = const; in fact, it is the value of J± when /3 ^ 7r/2 
which should be used in Eq. (0). In the pioneering study by SM, it was assumed that J± 
is initially zero. This assumption is valid only for grain temperatures approaching absolute 
zero. For nonzero grain temperatures the mean value of this component squared cannot 
be less than J]_q given by Eq. (^. To account for this non-zero component, while avoiding 
solving the corresponding Fokker-Planck equation, in our simplified treatment here we 
consider the evolution of the difference {J]_ — J±o), using the lucid approach suggested in 
SM. 
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In deriving Eq. (|38|) we assumed = const. Recognizing now that will be 
time-dependent, we note that disorientation of the grain depends primarily on the value of 
Jj^ near the time of crossover. We therefore seek to establish (J^(0))^/^, the value at the 
moment (t = 0) of crossover. One can write 

^i:^i-^ = iVi((AJ^r)-Mf^ , (40) 
at ts 

which generalizes eq. 37 in SM. All the time during a crossover, apart from a short interval 
when the grain actually flips over, u u!± and therefore it is possible to assume that ^ uj 
(see SM). According to our initial assumption, regular torques dominate the zero-order 
dynamics. Thus dt = {Iz/Lz) x = {Iz/Ni{AJz))duJz follows from Eq. (|36|) . Substituting 

c = (^^-) ( 2/. \ , _ , 



and 



/ 2/ \ 



into Eq. (0) gives (SM) 

1 - Cu'. (43) 



du 

Therefore for negative Uz increasing to zero for t = 0, one gets 



.3 

and 

,torque ' ±0 ' 

where 

. 3i/3r(4/3) UzJ.. 



C(0) = 3i/3r ( \) (44) 



3\i/3^ .4 



((AJx)|) + ^((AJ^)J) 



r(-)(A.ArO^/^|(AJ.)|-2/3/,^/3((AJ^)^)(l + x) . (46) 
and 

^ - N, ((AJx)S> • ' 



r(x) is the gamma function, and 
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is the value of uo-, such that the crossover time (IzUJz) / {Ni{AJz)) equals the relaxation 
time. During this time the number of torque events is \uJz\Iz/{AJz), and the product of 
this number over the mean squared increment of angular momentum per torque event, i.e. 
[((AJj_)^) + N2/Ni{{AJ±)'^)], provides the estimate for {J]_ — J]_q) in accordance with 
Eq. (HD. 



Comparison of Eq. ( ^Gj ) and eq. (42) in SM shows that the mean value of arising 
from stochastic torques is increased by a factor (1 + x), where 

and T is the gas temperature. In molecular gas with 1 — y <^ 1, x can be large, but in HI 
regions x ^ 0-5- We also note that x does not depend on grain geometry. 

Supersonic drift causes mechanical alignment that we briefly discuss in section 5.2. 
Here we limit discussion to the case where gaseous bombardment is isotropic during the 
crossover event. 

For typical interstellar conditions, the J]_q term in Eq. (^) is much more important 
than the term due to gaseous bombardment. The importance of the Barnett fluctuations 
relative to the stochastic torques is measured by 

7-2 m2/3 1 

= ^° ^ 412^ — ks{r) , (50) 

Jltorgue 7j/=^(l-|/)V3E2/32^1/6^'.j^l/3^5/3^1/3l + x 

^ 9^/3 X 4ir2 

M"") - (4^2_i)2/3[ig^2(^+i)+g^_^3] • (51) 



where 



The fact that the latter function tends to infinity for cubic grains (r = 1/2) is 
the consequence of the simplifications within our model. In fact, for cubic grains the 
perpendicular component of the grain angular moment will be of the order of the overall 
angular momentum, as pointed out in §2. It is easy to see, that for moderately oblate 
grains, however, {J]_{0)) is dominated by the term J^q arising from thermal fluctuations. As 
X ~ (1 — y)~^, for small concentration of atomic hydrogen J]_ torque become important. 
One of the problems with considering very small concentrations of molecular hydrogen is 
that disorientation then happens not only during crossovers but also during spin-ups (see 
Lazarian 1995d) and this requires the theory of suprathermal alignment to be modified. 
Moreover, according to Eq. ( [23| ) for y — >■ 1 only very large grains obey the theory we discuss 
here. 
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On estimating the critical size of in Eq. (^) we assumed that (J^) = J\q. In 
general, this is not true and our estimate of should be multiphed by (1 + R^'^Y^^ . The 
latter value, however, is ~ 1 according to Eq. ( [501 ) and therefore our estimate of the critical 
size given by Eq. (ESI) stays essentially unaltered. 



Although in the grain frame of reference J_L^torque and J±q appear very similar, their 
difference is obvious in the inertial frame. Indeed, J±o that arises from thermal fluctuations 
within the grain material does not alter the direction of J in the latter frame. On the 
contrary, J_L^torque that arises from gaseous bombardment and stochastic events of H2 
formation does directly affect the direction of J. 

As seen from Eq. (|50D , we expect to have J^q ^ J]_,torgueJ ^^i^ limit, the disorientation 
parameter F (see Eq. (|38|) ) can be obtained: 

F ^ 9.0 X W-^E'/^a'^^a-'/^f,'/'ps-'/\l + x)K,{r) , (52) 

where 

^ , /5 8(2r + 5)(4r^-l)^/' / 3(16r^(r + 1) + 6r + 3) ^ 

^^""^"VS 179ri/2(4r2 + 1)1/2 8r2(2r + 5) ) ' 

The disorientation decreases as r 1/2, which corresponds to a cubic grain. 
For the case of torque ^ ^lo "^^ "^^^ ^^^'^ obtain an estimate of F: 



F ^ 0.26 &^/^ps~^/^a~^/\hnE^^far^l^f^'^^^l + xK^{r) . (54) 

where 



^ 8 V205(2r + 5)(4r2- 1)1/6^1/2 / ^{iq^^^^ + I) + Qr + ?,)\ 

~ 179 3V6(1 + 4^2)i/2(ie^2(^ + 1) + 6r + 3)i/2 1^ + 8r2(2r + 5) ) ' 

(55) 

This estimate coincides with that in Lazarian (1995c) in the limit of negligible contribution 
from the gaseous bombardment^. 



4. Implications for the alignment 



^The difference in the numerical values obtained here and in Lazarian (1995c) stems from the fact that 
in the latter paper the function F was defined as the average value for an ensemble of grains with varying 
J_L and V (see §4.1). 
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4.1. Paramagnetic alignment 

Paramagnetic alignment of suprathermally rotating grains - frequently called Purcell 
alignment - can be described using the equation (Purcell 1979): 

d6 sin 6* cos 6* , ^, 

¥ = ' ^''^ 

where is the time of relaxation time of a grain with volume V in the ambient field B: 

where /sT ^ 1.2 x 10-^^77^ s (Draine 1996). 

The solution of the differential equation above is trivial: 

tan6' = tan6'oexp(— t/t^) . (58) 

If at t = grains are initially randomly oriented, then after time t we obtain the 
Purcell (1979) expression for 



^/ s , /r ,N 3 1 — (e'' — 1) arctan \/e^ — \ 1 , , 

g(t) = 3/2([cos2^]-l/3) = - ^ ^—3^ ^- , (59) 

where 8 = 2t/tr and square brackets denote averaging over grain initial orientations. 



Now suppose that grains are randomly oriented following crossovers and let P{t)dt be 
the probability that a randomly selected grain will have gone a time th G [t,t + dt] since its 
last crossover event. To obtain the Rayleigh reduction factor (Greenberg 1968) 

a = ^((cos2^)-l/3) , (60) 



one needs to average Q{t): 



S^Q{t)P{t)dt 



For our simplified treatment we will assume that for any particular grain in the 
ensemble the crossovers happen periodically with period tmax- Q Then 

P(Wx) = I J,"^^^; ^'"^^ ' (62) 



^Thc theory of alignment for an arbitrary distribution of time intervals between crossovers is given in 
Draine & Lazarian (1996). 
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For this distribution the mean time between crossovers (or zero-crossings) is = t 
Integrating ( |6TD we get 



max' 



max 



arctanV e'^™'"' — 1 ^ 



where 5max = ^tz/tj.. For small 5max EQ' (P^ can be expanded 



(63) 



5 52 53 

^ "^'^x "max I "max , (QA) 

10 210 840 ■ ^ ^ 

Up to now we assumed complete disorientation in the course of a crossover. It is evident 
from Table 2 that ^ for typical interstellar conditions. As the mean "time back to 
crossover" for "short-lived spin-up" (e.g. for < t^) is of the order of td, paramagnetic 
alignment is marginal unless the directions of J before and after crossovers are strongly 
correlated (SM). 

To account for incomplete disorientation SM adopted the following reasoning: consider 
crossovers that occur at intervals tmax! then in a time tmax/-^ the disorientation decreases 
cos?7 by 1/e. Thus, according to SM, the effects of incomplete disorientation during 
crossovers may be approximated by replacing by iz/mm[{F)j, 1]. Henceforth we use 
to denote averaging over the distribution of J±. We remind our reader that up to now we 

evaluated F for a grain with j± = {Jiy^^. 

We conjecture that replacing ^^ax in Eq. (pSf) 



(1 -exp(-(F)j)) 



to obtain 



3 

cr ^ 1 + — 

OeS 



arctanvt? 



Ve*<=ff - 1 

may give a better fit than the SM approximation above, as it allows for residual correlation 
for {F)j It is evident that for {F)j ^ 1 and {F)j <^ 1 our approximation coincides 
with that in SM. 

To study the effects of incomplete disorientation below we use Monte-Carlo simulations 
to calculate a for different ratios of iz/tr and {F)j. 

To obtain {F)j, we require (1/J^(0))j. To estimate this we note that 
{J]_{0)) = (^J(O)) + (^^(0)) and due to the symmetry inherent to the problem 



(4^(0)) = {J'ym = 0.5(Ji(0)) 



(67) 
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For a Gaussian distribution with af = (Jj(0)), one gets 



+~ 1 f rM , fn 1 a/tt 

-TT-I d^ = \/o- = 7727^ ■ (68) 



MO)/ J 27TafJo 2aj j V 2 ai (Ji(0))V 



Thus from Eq. (^) we get 



3/2 ((AJ)'^) 



4 ^KAJ.)KJi(0))V2 
= 1.60 X 10-2ii/2^i/2^-i/2^-i/2^^^-i/2(^ ^ X)i^5(r)(l + R-^)-^'\ (69) 

In those simulations we use Eq. (Q) to find {cos r]f) for a fixed J± and then perform 
numerical averaging of {cos rjf) over a Gaussian distribution of J±. We require a distribution 
function for the stochastic jumps that correspond to ((cos?7/))j. We assume the distribution 
of 7] to have the form^ 

Pi{ri)dri = C smr]exp{—a'^ sm'^{ri/2))dr] , (70) 
where C is the normalization constant: 



1 

Tl-exp(-a2^ ^'^^^ 



and a is the solution of the transcendental equation: 

1 + exp(-a2) 2 

(cos?7/ = . 72 

1 — exp(— a'') 

An individual jump over rj during a crossover event happens in a random direction and 
we obtain the final value of 6ij after the i-th crossover from the following formulae: 

cos 9ij = cos Oi^h cos T) + sin Oi^b sin rj cos x , (73) 

where x is a random variable uniformly distributed over [0, 2tt] and 6i^b is the value of the 
alignment angle just before the i-th crossover. Between crossovers the dynamics of the 
alignment angle is determined by Eq. (^S|). Averaging over P(tmax) (see Eq. (p^)) is also 
performed to account for the distributions of the times since the last crossover. 



^ The assumed functional form ( [70| ) has the required behavior of dPi/drj = for 77 = 0,7r; Pi ~ (sin 77) 
for a ^ 0; Pi ^ 77 exp(— q;^?7^/4) for r/ <C 1. 
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The results of those calculations are shown in Fig. 2 where we have plotted a vs S^.^ 
(defined by Eq. (p5|)). For each value of (5efr different symbols correspond to the alignment 
measures obtained for different values of {F)j. The solid line in the same plot corresponds 
to Eq. (^). In the limit {F)j —^oowe have complete disorientation, in which case Eq. ( |5BD 
is an exact result (for periodic crossovers). However, it is evident that Eq. (|66|) provides a 
good approximation to the numerical results for finite (-F)j, at least for periodic crossovers. 
More general models where the times between crossovers are obtained through Monte-Carlo 
simulations are studied in Draine & Lazarian (1997). 

For typical values of interstellar parameters (see Table 2) one obtains F 0.014 (see 
Eq. ^2p. Using our earlier estimate ~ LGt^d^/^ (for assumed ti/td = 0.25/a) we get 
6es ~ 13.7 (indicated by an arrow in Fig. 2), for which Eq. (|66| ) gives a ~ 0.8. This high 
degree of alignment is due to the small estimated value of (F). In fact, it was argued in 
Lazarian (1995c) that is expected to be several times greater than for grains with 
a > 10~^ cm and this, by increasing tz, would further increase the expected alignment. 

Although we do not quantitatively discuss here the alignment of grains with 
a < Oc ~ 1.5 X 10~^ cm we conjecture that the alignment of such grains may be suppressed 
by the effective disalignment during each crossover (i.e. (F) ^ 1), due to the variations in 
the angle j3 due to Barnett fluctuations when <^ (see §3.2). The strong dependence of 
ts/tc on a (see Eq.(p2D) suggests that this may account for the observed lack of alignment 
of interstellar grains with a ^ 10"^ (Kim & Martin 1995). 

We see, then, that if the only important torques were those due to H2 formation, 
gas-grain collisions, and paramagnetic dissipation, we would expect paramagnetic grains in 
diffuse clouds to be substantially aligned for a > 0^ ~ 1.5 x 10~^ cm, and probably minimally 
aligned for a < ac, at least qualitatively consistent with observations. It has recently been 
recognized, however, that starlight plays a major role in the dynamics of a 0.1 /im grains: 
the torques exerted by anisotropic starlight (z) drive suprathermal rotation (Draine & 
Weingartner 1996) and (u) can directly act to align J with the interstellar magnetic fleld 
(Draine & Weingartner 1997). 

Since we have neglected starlight torques in the present paper, our conclusion for 
a 0.1 /xm grains is only preliminary. A study of crossovers incorporating the effects of 
both H2 formation and anisotropic starlight is planned. It appears to us highly likely that 
when both effects are included, the observed lack of alignment of a ^ 0.1 /im grains, and 
substantial alignment for a 0.1 fim grains, will be explained. 

We recall that suprathermal rotation can also be driven by variations of the 
accommodation coefficient and photoelectric emissivity (Purcell 1979). In a molecular 
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{y = 1) region with no ultraviolet light, Purcell's estimate for the torque due to variations 
in accommodation coefficient leads to Oc = 3.0 x 10^^ cm as the radius for which ts = tc- 
However, for this case we also find (F) > 1 for a ~ Oc, so that Purcell alignment will be 
insufficient unless iz > tr- 

4.2. Mechanical alignment 

It was previously thought that suprathermally rotating grains are not subject to 
mechanical alignment when the gas is streaming relative to the grain. However, two 
mechanisms of mechanical alignment of suprathermally rotating grains, namely, "cross- 
section" and "crossover" alignment, were proposed by Lazarian (1995d). The first process 
is caused by the fact that the frequency of crossover events depends on the value of the 
cross-section exposed to the gaseous fiux (see Lazarian & Efroimsky 1996 for more details). 
The second mechanism arises from the substantial susceptibility of grains to alignment by 
gas-grain streaming during crossover events. 

Both mechanical processes are related to the phenomenon of crossovers. Thus 
our finding of reduced disorientation during crossovers is a new feature that should be 
incorporated into the discussion of mechanical alignment. As we mentioned earlier, this 
reduced randomization is valid only for grains with a > a^, where Oc is given by Eq. (|23|) , 
and therefore no changes of the earlier results are expected for grains with a < 10~^ cm. 
Such grains can be aligned, for instance, by ambipolar diffusion, which favors small grains. 

To start with, consider the cross-section mechanism. In Lazarian (1995d) this 
mechanism was exemplified using a toy model, namely, a fiat disc grain which randomly 
jumps in the course of a crossover from one position, where the surface of the disc is parallel 
to the fiow, to the other position, where the disc surface is perpendicular to the fiow. If the 
probability per unit time of a crossover is proportional to the rate at which atoms arrive at 
the grain surface, it is easy to see that the grain will spend more time at orientations where 
the cross section presented to the streaming gas is minimal. Within the toy model above, 
this corresponds to the position with the surface of the disc parallel to the fiow. 

In other words, the cross-section mechanism uses the fact that is a function of the 
angle between the grain axis of major inertia and the direction of the gaseous fiow. 
Roughly speaking, our study above shows that crossovers with disorientation parameter 
(F) and mean time between crossovers tz are equivalent to crossovers with complete 
disorientation and the mean time between crossovers iz/{l — exp(— (F))). If F is dominated 
by thermal fiuctuations its dependence on gaseous bombardment vanishes (see Eq. (^)) 



- 20 - 



as does its dependence on 0. Thus the only effect of incomplete disorientation during 
crossovers (as compared to full disorientation) is to increase of alignment time (the time to 
attain a steady-state) by a factor (1 — exp(— (F)))~^. 

"Crossover alignment" depends on the ratio of the randomizing torques arising from H2 
formation and aligning torques caused by gaseous bombardment (see Lazarian 1995d). This 
ratio neither depends on the number of crossovers nor on the time of alignment. The fact 
that the thermal fluctuations do not change the direction of J is essential for understanding 
why this type of alignment is not suppressed in the presence of the incomplete disorientation 
during crossovers. It is possible to show, however, that the time of alignment increases by a 
factor (1 -exp(-(F)))-^ 

In spite of the fact that the measure of the mechanical alignment does not change, 
our observation that the time required to reach steady state is increased by a factor 
(1 — exp(— (F)))^^ can be important. This is particularly important whenever grain 
alignment is caused by a transient phenomenon, e.g., a MHD shock. If ^^/(l — exp(— (F))) 
is much longer than the time of streaming, the alignment of grains with a > 1.5 x 10~^ cm 
will be marginal^. 

5. Conclusion 

We have shown that thermal fluctuations within the grain material limit the extent 
to which the axis ai of major inertia can be aligned with the angular momentum J in 
suprathermally rotating grains. Although the fluctuating angle (3 between ai and J is tiny 
when the grain is rotating suprathermally, it becomes larger and of critical importance 
during periods of crossovers. We have proved that for grains with a > Oc ~ 1-5 x 10"^ 
cm the non-zero component of ai perpendicular to J arising from thermal fluctuations 
substantially diminishes the degree of the randomization of the angular momentum direction 
in the course of crossovers. If the only torques acting on a grain are those due to gas-grain 
collisions, H2 formation, and paramagnetic dissipation, our estimates show that for large 
(a ^ cic) grains the grain alignment is close to perfect, while small (a < Qc) would have only 
marginal alignment. 

If there is gas-grain streaming, the thermal fluctuations increase the time for mechanical 
alignment for large suprathermally rotating grains, but do not alter the limiting steady 



^This by no means preclude small grains from being aligned and we believe that the preferential alignment 
of small grains could be a signature of such an alignment. 
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state measure of alignment. If the mechanical alignment is caused by Alfvenic waves, it acts 
in unison with the paramagnetic mechanism to enhance the alignment of large grains. For 
small grains mechanical alignment due to transient phenomena (e.g. ambipolar diffusion 
within MHD shocks) can be the dominant cause of alignment. 

A.L. is extremely grateful to Lyman Spitzer for elucidating discussions of the crossover 
effect. A.L. acknowledges the support of NASA grant NAGS 2858 and B.T.D. the support 
of NSF grant AST-9219283. 



A. Some Results for a Square Prism 

Here we consider a square prism, with dimensions b x b x a, density ps, 
mass Psb'^CL-i area 46a + 26^, and moments of inertia = {8/3)psa{b/2)^ and 

= = (l/3)p,a(6/2)2[a2 + 4(6/2)2]. We let r = (6/2a) and note that r = 1 
corresponds to the prism grain discussed in Purcell (1979). We consider a hydrogen-helium 
gas with density nn = n(H) + 2n{E.2)i temperature T, molecular fraction 

; (Ai) 

and ?T,(He)/?T,H = 0.1. The square prism has 

pstt 2r (r + 1) 

™mH(1.2 -0.293y) (r + 1) ' ^'=(7+2)^*' ' ^ ' 

where tM is the time for the grain to collide with its own mass of gas, td is the rotational 
damping time (assuming incident atoms to temporarily stick) ,[3 and vu = {SkT/nmuY^'^ is 
the mean speed of H atoms. If a fraction 7 of impinging H atoms are converted to H2, then 
the H2 formation rate is 

Ni = r(r + 1)7(1 - y)nnvua^ • (A3) 

We assume the grain to be spinning around the 2;-axis. The prism is assumed to have u 
active sites of H2 formation distributed randomly over the surface. Following Purcell, we 
assume that newly-formed H2 molecules depart from each recombination site at a rate 
Ni/u, with fixed kinetic energy E but random directions [dP/dO = 2 sin ^ cos ^, where 6 is 



10 pm-(2gii ^ Spitzer (1971) give tj_/tM for a square prism but their eq.(9) contains a typographical error: 
the factor (5s + 1) should instead be (4s + 1). Our r is equal to l/(2s) as defined in Purcell & Spitzer (1971). 
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with respect to the local surface normal). The v/{l-\-r) sites on the sides of the prism then 
produce a steady torque with 

= + 1)^/^(1 -y)nH^Ha^(^)'^' , (A4) 
and a mean angular impulse per recombination event (AJ^), with 

KAJ.)|~ - + 1,1/2 3^ J ■ 

Individual H recombination events, occurring at a rate Ni, contribute random angular 
momentum impulses with 

{{/^mn = ^^^^y^rnWE , (A6) 

u^J ^2^ 16r2(r + l) + 6r + 3 , 
((AJ^)^)h = i2{r+l) ""^"^^ ■ ^^^^ 

Gas particles impinge at a rate 

A^a = 2r(r + l)nHVHa2(1.05-t/ + y/V8) , (A8) 

If impinging particles temporarily stick and then thermally desorb at temperature 7^, then 
these collision events produce 

N2{{AJ,)^)g = '^r\2r + h)nnm^vnxa%T + Ta) (l.2 - ?/ + ^) > (A9) 

N2{{AJi_f)g = U{lQ{r\r + 1) + 6r + ^)n^mY,v^a''k{T + T^) (^1.2 - ?/ + ^) • (AlO) 
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notation 


meaning 


y = 2n{H2)/nH 


H2 fraction 


7 = 0.27 


H recombination efficiency 


vh = {SkTgas/T^mHf^ = 1.5 X lO^VH cm s"^ 


thermal velocity 


a = 10~^ cm 


grain size 


r = b/2a 


grain axis ratio divided by 2 


E = 0.2E eV 


kinetic energy of nascent H2 


a = lO^^a cm~^ 


surface density of recombination sites 


V = 80r(r + l)ad'^ 


number of recombination sites 


Ta = Ibfd K 


grain temperature 


T = 85fK 


gas temperature 


Uh = 20nH cm"^ 


density of H nucleon 


Qs = ^Qs g cm-3 


solid density 


B^bB X 10-^ G 


magnetic field 



Table 1: Parameters of grains and ambient medium adopted in this paper. 



Crossover = 2.9 x lO^ipTduE-^d-^y/^ /{^nnVHil - y))Zo{r) s 

Barnett effect ts = 7.1 x lO'^d'^ p1{uj_io/u;yZi{r) s 

Gaseous damping — 3.3 x 10^'^dps{n?HifH)^^{l.2 — y + y/\/2)Z2{r) s 

Paramagnetic relaxation tr = Q.7 x 10^^ Psd'^TdB~'^Z3{r) s 



Table 2: Characteristic times involved. We use K 1.2 x 10 ^^T^ ^ for paramagnetic 

relaxation (Drainc 1996). The functions of grain axis ratio Zi i = 0,3 arc as follows: 
Zo(r) = (6(l+4r2))V2/(5(4r2-l)(r + l))V2, Zi{r) = {4r^ + if/ {25r^) , Z^ir) ^ 2r / {r + 2) 
and Z^{r) = r^. 




Fig. 1. — Grain's axis of major inertia z 
momentum J. 



and precess about the direction of angular 




Fig. 2. — Measure of alignment a [see Eq. (0)] for periodic crossovers as a function of Jeff, 
defined by Eq. (pS]). The solid line is the analytic estimate (p^).The arrow corresponds to 
(5eff estimated for typical interstellar conditions (see text). 



